Confinement resonances in photoionization of endohedral atoms: a myth or reality? 
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We demonstrate that the structure of confinement resonances in the photoionization cross section 
of an endohedral atom is very sensitive to the mean displacement (a) of the atom from the cage 
center. The resonances are strongly suppressed if 2(a) exceeds the photoelectron half- wavelength. 
We explain the results of recent experiments which contradict the earlier theoretical predictions on 
the existence of confinement resonances in particular endohedral systems. 
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In this Letter we explain the discrepancy between the- 
oretical predictions and experimental results for the pho- 
toionization (PI) of endohedral atoms A@Cjv- The dis- 
crepancy concerns 'confinement resonances' [I] in the PI 
spectrum whose existence was predicted theoretically [2- 
17] (and references therein) but so far has been not sup- 
ported experimentally [18-20]. We also formulate the 
criterion which allows one to estimate, for a particular 
endohedral system, the interval of photon energies where 
the resonances can appear. 

Confinement resonances appear as a result of interfer- 
ence between a direct wave of the photoelectron escaping 
the atom and the waves due to scattering from the atoms 
of the cage. Depending on the photoelectron momen- 
tum the interference can be constructive or destructive. 
Thus, the spectrum of the encaged atom acquires addi- 
tional oscillations as compared to the free atom. It was 
noted [2] that the cage- induced oscillations have the same 
nature as the extended X-ray absorption fine-structure 
(EXAFS) for solid-state systems. 

To date, a number of theoretical investigations have 
been carried out of the features of confinement reso- 
nances for various endohedral systems. The resonances 
were studied for noble-gas [2, 5, 7, 12, 15] and metal [2- 
6, 8, 9, 16] atoms encaged in Ceo, for Ar encaged in larger 
fullerenes C240 and C540 as well as in onion-like structures 
C6o@C24o etc [15], for hollow atoms [14]. A strong dis- 
tortion of the atomic 4d-f giant dipole resonance due to 
the scattering from the cage was predicted for Ba@C6o 
[2, 3], La@C 60 [2], Xe@C 60 [12] and Ce 3+ @C 60 [16]. The 
confinement resonances were also investigated for the an- 
gular distribution of photoelectrons [6, 10, 11]. Specific 
features of the resonances in the PI of endohedral ions 
A@C^ 2 were studied in Ref. 13. These and other closely 
related topics can be found in the recent review [17]. 

The experimental data on the PI of A@Cjv are sparse 
due to '. . . the difficulty to produce sufficient amounts of 
purified endohedrals for the gas phase experiments' [20]. 
So far the PI cross section in the region of the 4d-f giant 
resonance were measured for Ce@Cs2 [18], Pr@Cs2 [19] 
and Ce@C^" 2 [20]. In [20] it is explicitly stated that no 
confinement resonances were detected. With respect to 



[18, 19] one draws this conclusion matching the measured 
data to big absolute uncertainties. 

Most of the theoretical work has been devoted to the 
confinement resonances in spherical endohedral systems, 
in which the atom is placed at the center of a spheri- 
cal cage. The correlation between the confinement res- 
onances and the off-the-center position of the atom was 
reported for Ba@C 60 [4], Mg@C 60 and Mg@C^ [9]. A 
strong sensitivity of the photoelectron angular distribu- 
tion and total cross section of PI on the displacement 
of the atom from cage center was demonstrated [8, 10]. 
Specific results were presented for Li@Ceo and Ar@C6o- 

In the cited papers the influence of the off-the-center 
position was analyzed for the fixed-in-the-space endohe- 
dral fullerenes, i.e. fixing both the direction and magni- 
tude of the displacement vector a between the atom and 
the center. In experimental conditions cither one of these 
parameters or both (depending on the endohedral sys- 
tem) cannot be controlled. As a result, proper averaging 
procedures must be adopted in order to bring theoretical 
predictions closer to the observable quantities. 

For example, in metallofullerenes M@Cjv (M stands 
for a metal atom), the atom, as a rule, is located off the 
center (see, e.g., Ref. 21). The equilibrium can be very 
stable due to a strong bond of M to the carbon atoms. 
Hence, one can fix the distance a = |a|. To explore the 
properties of an individual metallofullerene the direction 
of a also can be fixed. However, to compare with the gas 
phase experiments, any observable, including the cross 
section, must be averaged over the directions. In this 
connection we mention Ref. [6] where the atomic Is PI 
of M@C6o molecules (M=Li, Na, K) was theoretically 
investigated for the gas of fullerenes. The authors, in- 
dicating the off-the-center position of the atom, ignore 
its influence on the cross section and carry out the study 
placing the atom at the center for 'computational sim- 
plicity' (see the end of Introduction). In what follows 
we demonstrate, that, generally, such an assumption is 
not correct and can lead to erroneous predictions on the 
existence of the confinement resonances. 

Other complexes of current interest are endohedral 
noble-gas atoms [17]. In these systems the atom stays 
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neutral, and its dynamics is defined by the Lennard-Jones 
interaction with the cage atoms [22, 23]. The equilibrium 
position and the amplitude of thermal vibrations depend 
on the atom type, the cage size and the temperature. 
For such complexes it is reasonable to carry out also the 
averaging over the a values. 

As we demonstrate below, in many cases the averag- 
ing destroys the confinement resonance structure. Prior 
to introducing the formalism, let us present qualitative 
arguments. For simplicity, we ignore multiple scattering 
from the cage. Then, the oscillations in the cross section 
arc due to the interference of the two waves of photoelec- 
tron, the direct and the scattered ones, both originating 
from the same source, - the encaged atom. The atom at 
the center can be treated as a point-like source. For a / 
the source acquires a size D ss 2(a) where (a) is the mean 
distance from the center. The finite size of the source in- 
fluences the interference pattern. When D exceeds the 
half-wavelength the pattern is destroyed. Hence, one can 
expect that for the photoelectron momenta [24] 
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the confinement resonances disappear. The correspond- 
ing range of photon energies is uj > Lu m - U1 — pf nin /2 + Iq, 
where Iq is the ionization potential of the atomic shell. 

To start with the formalism let us note that, since the 
confinement resonances are due to the interference phe- 
nomenon, one can adopt the simplest model for electron- 
fullerene interaction which ensures the interference. We 
model this interaction with a <5-like potential well [5, 25]. 
Choosing the origin at the atom and introducing the po- 
sition vector a of the cage center, one writes: 



C/ S (r;a) = -A5(\t-sl\ - R) , 



(2) 



where A can be related to the cage radius R and the elec- 
tron affinity of the fullerene [5, 25]. For the central posi- 
tion of the atom eq. (2) defines the spherically-symmetric 
potential U s (r) = —AS(r — R) used in the cited papers. 

The wavefunction \fp ^ (r; a) of the photoelectron emit- 
ted from the encaged atom satisfies the equation [26]: 

¥-\r; a) =^ p \r) -jdv'G £ {v, r')U s (r'; a)* p -)(r'; a).(3) 

The superscript (— ) indicates that the asymptotic form 
is 'plane wave + incoming wave', p and £ = p 2 /2 are 
the momentum and energy, ijj p ^(r) and G e (r, r') are the 
wavefunction and the retarded Green function of the elec- 
tron escaping from the free atom. 

For an at-the-center atom eq. (3) is solved by expand- 
ing "fp '(r; 0), tpp \r) and G £ (r, r') in series over spheri- 
cal harmonics. Then, evaluating the integral and match- 
ing both sides of the equation at r = R, one derives 



where a = 27rA, and Yj m are spherical harmonics. The 
radial wavefunctions P p i(r) and the scattering phaseshifts 
Si = Si (p) correspond to the motion in the field of atomic 
residue. The radial Green function g p i{r,r') is propor- 
tional to the product of the regular P p i (r) and irregular 
X P i (r) (at r = 0) solutions of the radial Schrodinger equa- 
tion: g p i{r,r') = -P pi (r<) x P i(r>)/2p with r>/r< being 
the largest /smallest of r,r' (see, e.g., [26]). 

Using (4) one relates the amplitude and cross sec- 
tion of PI of the cndohcdral atom to those of the free 
atom. For the free atom, the partial amplitude M A 
of the dipole transition from the subshell vq = (jiqIq) 
(jiq is the principal quantum number, Iq is the orbital 
momentum) is proportional to the matrix element of 
the radial dipole operator d: M^ q oc J q °° P p i (r) d P Vo (r) 
where P» a (r) is the shell wavefunction and I = Iq ± 1 
due to the dipole selection rules. To evaluate the am- 
plitude for the encaged atom one substitutes P p i (r) with 
the expression in the brackets in (4). Additionally, as- 
suming the shell radius to be much smaller than R, 
one sets P p i (R) g p i (r, R) = P p i (r) g p i (R, R) and derives 
M A@Cjv (a = 0) w M;f o /(l + ag p i{R,R)). The ratio of 
the partial cross sections cr^® CjY (w; a) / aft (u>) = ?// (oj; a), 
calculated for a = 0, reads 
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■g p i(R, R) 



l-2aReg p i(R, R) . (5) 



The photon energy lu is related to p via p = \J2(oj + Iq). 
The first relation in (5) is similar to the formula derived 
in [5, 25]. The approximate relation is valid when the 
cage potential is treated perturbatively. 

For momenta p large enough that pR 3> 1, one can use 
the asymptotic formula g p i(R,R) oc cxp(— ipR) sin(pi? — 
ttI/2 + 5i). Then, eq. (5) explicitly reveals the oscillatory 
behavior of t]i(uj; 0) as a function of p (or ui) due to the 
interference of the direct and scattered waves. 

The exact solution of eq. (3) with the potential (2) 
taken for the off '-the- center position has not been found. 
However, as mentioned, the cage potential can be treated 
perturbatively. The first-order solution of (3) reads 
* p _) (r;a) w ipj~ \r) + Aifj ( p ~\r; a) , where 



AV> p _) (r,a) =A/dr'G(- ) (r,r')5(|r'-a| - R)tp. 



p- } (r')(6) 



xy lm (n r )y ( ;(fip), 



(4) 



is the scattered wave. The total amplitude of PI is 
calculated as M a@Cn (a) = /dr* p _) *(r; a) e • d>o( r ) = 
M A + AM{a), where M A is the amplitude for the free 
atom and AM(a) is due to (6). Other notations include 
the unit vector of the photon polarization e, the oper- 
ator of the dipole moment d, the wavefunction of the 
subshell tpo(r). Squaring the modulus of the amplitude 
and retaining the terms linear in the potential, one writes 
the cross sections as ct a@Cjv (w; a) a A (oj) + Aer(u;;a), 
where the term Aa(u>;&) oc Rc(.M A AX*(a)) contains 
the dependence on a and, thus, must be averaged. 
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First wc consider the angular averaging cr A@Cjv (u; a) = 
f <t a@Cn (uj;&)dn a /4:iT. The dependence on the direc- 
tion of a enters A<7(a;;a) via the 5-function from (6). 
For its averaging one can use the expansion (e.g., [27]): 
5(|r - a| - R) = (2nR/ar) £ im Pi(0^(^)^(n a ), 
where Pj(£) are the Legendre polynomials, and £ = (r 2 + 
a 2 - R 2 )/2ar subject to |£| < 1. Hence <5(|r - a| - R) = 
R/2ar if £| < 1 and = if otherwise. The rest of the alge- 
bra is straightforward. To evaluate the radial integral in 
AM (a) one uses the approximation discussed in connec- 
tion with eq. (5). Finally, noticing that for each partial 
transition Iq — > I = Iq ± 1 the averaged term A.a(u>; a) is 
one derives the following expres- 

A 



proportional to ajj (u 



sions for the ratio r] l (ui]a) = o 11q n (uj\ a)/a^ Q (uj 



m{u;a) = 1 - a 



R+a 

R f d ^ Reg p ;(r, r) 



(7) 



R-a 



1 + (-1)'q S p {a)cos2Si + C p {a)sm25i 



(R/2pa) Jr^q dr r 1 sin(2pr) and 



where S p (a) — ^^~ y u,j jR _ a 

C p (a) = (R/2pa) drr' 1 cos(2pr). 

The first relation in (7) qualitatively explains the im- 
pact of the angular averaging. For a = the formula 
reproduces the right-hand side of eq. (5), derived for the 
centrally positioned atom treated as a point-like source. 
As a increases the averaging leads to a non-zero effec- 
tive size of the source, D = 2a. When D becomes larger 
than the emitted half- wavelength the interference is lost. 
Hence, for p > p m i n = ir/2a the confinement resonances 
in the cross section profile disappear. 

The right-hand side of (7) allows one to analyze quanti- 
tatively the modification of the interference pattern with- 
out calculating the parameters of photoclcctron wave- 
function (in particular, the phaseshifts) but using only p, 
R and a. Indeed, the relative change in the amplitude 
of the oscillations can be understood by comparing the 
functions S p (a) and C p {a) to their values at a = 0. 

Let us apply eq. (7) to estimate the effect for Ce@Cs2- 
In [18, 20] the PI cross section was measured in the region 
of the 4d giant dipole resonance, i.e. for w = 120 . . . 140 
eV, and no oscillations were seen. The 4d ionization po- 
tential in Ce is 114 eV for 4d3/2 shell and 111 eV for 4ds/2 
shell [28]. Hence, the indicated range of to corresponds 
the photoclcctron momenta p = 0.5 ... 1.5 a.u. The mean 
radius of Cs2 is 4.15 A and Ce atom is displaced by 1.8 A 
from the center [29]. The dependences S p (a) and C p (a) 
on p arc presented in Fig. 1. It is clearly seen that angular 
averaging destroys the interference for p > p m i n ~ 0.46 
a.u, i.e. for the photon energies u > Iq + 3 eV. 

Similar estimates explain the absence of the confine- 
ment resonances for Pr@Cs2 [19]. Less univocal conclu- 
sion can be drawn for Li@Cgo, K@C6o and Mg@C6o the 
PI of which was considered theoretically [6]. The au- 
thors indicate that the atoms displacements arc 1.5, 0.25 
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FIG. 1: S p (a) and C p (a) versus p calculated for Ce@Cg2 for 
at-the-center position of Ce and off-the-center with a= 1.8 A. 



and 1 A for Li, K, and Mg, respectively, but set a = 
when calculating the cross section. As a result, for each 
metallofullerene they predicted several confinement reso- 
nances for the photoelectron energies below 2 a.u. How- 
ever, calculating e m i n = p 2 nin /2 and obtaining the values 
0.15, 5.5 and 0.35 a.u., we state that the resonances in 
Li@Cgo and MgQCgo will hardly survive the averaging 
procedure, whereas the resonances predicted for K@Cgo 
will be preserved due to small value of a. 

To develop the theory further one accounts for ther- 
mal vibrations of the atom in the vicinity of equilibrium 
position. This can be important, e.g., for endohcdral 
noble-gas atoms. Let us construct a model potential for 
the A-Cn interaction which, despite being quite crude, 
allows one to analyze the influence of the vibrations on 
the interference pattern. Ignoring the atom-fullcrcne hy- 
bridization one builds the potential as a sum of pairwise 
A-C Lcnnard- Jones potentials over all carbon atoms [22]. 
Assuming a homogeneous distribution of carbon atoms 
over the sphere of radius R and substituting the sum 
with the surface integral, one derives the A-Cjv potential 
U (a) = Ui,(a) + U CYi (a) as a function of a [30]. The term 
Ui^ia) is the attractive (the London-type) potential, and 
U ex (a) is the repulsive exchange term. 
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FIG. 2: Potentials U(a) for various A@Cjv as indicated. 



The functions U (a) for various systems are presented 
in Fig. 2. The radius 3.5 A of the Ceo cage is less (for Ar, 
Kr, Xe) or just above (He, Ne) the equilibrium distance 
of the Lennard- Jones potential [22], so that the atom 
is confined to the center [31]. For C240 the radius 7.1 
A exceeds the equilibrium distance. Hence, the atom is 
either weakly bound well off the center (Ar, Kr, Xe) or 
moves freely in most part of the fullerene (He, Ne). 

The probability to find an atom at the distance a is 
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given by dWr(a) = Cexp (— U{a)/kT) a 2 da, where T is 
the temperature, /c - the Boltzmann constant, and C = 
[/ a * nax dWr(a)] with a max ~ R — Ra {Ra is the radius 
of the atom). To average the ratio rji(uj;a) over a one 
multiplies eq. (7) by dWr(a) and integrates over a: 

foH) = 1 + a {(S p ) cos 2(5/ + (C p ) sin 2,5,] 

with the general notation {g) = Jq™"* x g(a) dWr(a). 
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FIG. 3: versus p calculated for various T and A@Cjv. 

Dependences {S p ) on p are presented in Fig. 3 [32]. It 
is seen that in all A@C24o systems the oscillations are 
suppressed because of strong off-the-center position of 
the atom. The situation is not so universal for a smaller 
Ceo fullercnc in which the equilibrium is at the center. 
Indeed, for Ne the amplitude of oscillations decreases no- 
ticeably for all T starting with p « 1 a.u., whereas for 
Xe the oscillations stay nearly unchanged up to p > 1 
a.u. To explain the difference one estimates the effective 
size of the source as the doubled mean amplitude of ther- 
mal vibrations and calculates the momenta p m i n = n/D 
beyond which the interference disappears. For T = 500 
K the calculated values of p m ; n are 1.8, 3.9, 7.1 a.u. for 
Ne@C 60 , Ar@C 60 and Xe@C 60 , respectively. 

In summary, we have demonstrated that the interfer- 
ence effects in photoionization of A@Cjv arc very sen- 
sitive to the mean displacement (a) of the atom from 
the cage center. The range of photon energies, where 
the confinement resonances can be observed, one deduces 
by matching 2(a) the photoelectron half- wavelength. In 
metallofullerenes, due to a non-central position of the 
atom, the angular averaging can destroy the oscillatory 
structure predicted for the atom at-the-center. For large 
noble-gas endohedral fullerenes (C240 and larger) the res- 
onances do not survive the averaging procedure with 
the Boltzmann distribution function. For a smaller Ceo 
fullerene a more rigorous treatment of the atom- fullercnc 
interaction is needed to draw the final conclusion on the 
existence of the resonances. The role of non-centrality 



must be studies for other related phenomena (Coulomb 
confinement resonances, non-dipolc effects, the photoion- 
ization of onion- like structures etc [17]). 
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